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Abstract 

Nonlinear boundary value problems (BVPs) by means of the classical Lie symmetry 
■ method are studied. A new definition of Lie invariance for BVPs is proposed by the 

generalization of existing those on much wider class of BVPs. A class of two-dimensional 
nonlinear boundary value problems, modeling the process of melting and evaporation of 
_ metals, is studied in details. Using the definition proposed, all possible Lie symmetries and 

Q ■ the relevant reductions (with physical meaning) to BVPs for ordinary differential equa- 

. tions are constructed. An example how to construct exact solution of the problem with 

correctly-specified coefficients is presented and compared with the results of numerical 
^Njj' ■ simulations published earlier. 

2010 Mathematics Subject Classification : 22E70, 35K61, 80A22. 

-S' 1 Introduction 

It is well known that principle of linear superposition cannot be applied to generate new 
exact solutions to non/mear partial differential equations (PDEs). Thus, the classical methods 
CN I (the Fourier method, the method of the Laplace transformations, and so forth) are not applica- 

\^ • ble for solving nonlinear PDEs. While there is no existing general theory for integrating non- 

O ■ linear PDEs, construction of particular exact solutions for these equations is a non-trivial and 

important problem. Now the most popular methods for construction of exact solutions to non- 
integrable nonlinear PDEs are the Lie, Lie-Backlund and conditional symmetry methods [IHS] . 
Although these methods are very powerful provided the relevant symmetry is known, several 
other approaches for solving non-integrable nonlinear PDEs were independently suggested dur- 
ing the last decades. Among them the method of compatible differential constraints [HIH], the 
method of linear invariant subspaces [9], tanh- method and its various modifications [TOl - ITi] . 
^ ! the method of additional generating conditions [15l|T6] , and the transformed rational function 

c3 i method [T7] should be marked out (see, e.g.. Supplements in [IB] about other methods). 

The Lie symmetries are widely applied to study nonlinear differential equations (including 
multi- component systems of PDEs) since 60-s of the last century, notably, for constructing their 
exact solutions. Nevertheless there are a huge number of papers and many excellent books (see, 
e.g., P-|5 and papers cited therein) devoted to such applications, one may note that a very 
small number of them involve Lie symmetries to solve boundary value problems for the given 
PDEs. To the best of our knowledge, the first papers in this directions were published in 
the beginning of 1970-s [12] and [20] (the extended versions of these papers are presented in 
books pi] and [2], respectively). The books, which highlight essential role of Lie symmetries 
in solving boundary value problems (BVPs) and present several examples, were published only 
in 1989 [2l[22]. 

The main object of this paper is a nonlinear EVP of Stefan type, which belongs to the class 
of BVPs with free (moving) boundaries. Boundary value problems of Stefan type are widely 
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used in mathematical modeling a huge number of processes, which arise in physics, biology 
and industry [231 - I27] . Nevertheless these processes can be very different from formal point of 
view, they have the common peculiarity, unknown moving boundaries. Movement of unknown 
boundaries are described by famous Stefan boundary conditions [271I2S]- It is well-known that 
exact solutions of BVPs of Stefan type can be derived only in exceptional cases and the relevant 
list is rather short at the present time (see [23 t [291435] and papers cited therein). 

Nevertheless BVPs with free boundaries are more complicated objects than the standard 
BVPs with the fixed boundaries, it can be noted that the Lie symmetry method should be 
more applicable just for solving problems with moving boundaries. In fact, the structure of 
such boundaries may depend on invariant variable (s) and this gives a possibihty to reduce the 
given BVP to that of lover dimensionality. This is the reason why different authors applied 
the Lie symmetry method to BVPs with free boundaries ignoring BVPs with fixed boundaries 

[IHlEolESHsH]. 

The paper is organized as follows. In Section 2, we discuss the existing definitions of 
Lie invariance for BVPs and propose their generalization on much wider class of BVPs. As 
an example the direct application of the definition for the well-known BVP with the fixed 
boundaries is presented. In Section 3, we apply the definition derived to the class of (1+1)- 
dimensional BVPs of Stefan type used to describe melting and evaporation of materials in the 
case when their surface is exposed to a powerful flux of energy [311 [39]. The result obtained 
is an essential generalization of paper [33 . In Section 4, we reduce the problem to BVPs for 
ordinary differential equations, using Lie symmetry operators obtained in the previous section. 
An example how to construct exact solution of the problem with correctly-specified coefficients 
is also presented. Finally, we present conclusions in the last section. 



2 Definition of Lie invariance for BVPs 

We start from a definition of invariance of a BVP under the given infinitesimal operator 
presented in [21 [5] and restrict ourselves on the case when the basic equation of BVP is an 
(l+l)-dimensional evolution PDE of /cth-order {k > 2). In this case the relevant BVP may be 
formulated as follows: 

Ut = F (x, u,u^,..., u^^^) , (t, x) G ^] C (1) 

Sa{t,x) = 0: Ba{t,x,u,u^,...,u^^~^^) =0, a=l,2,...,p, (2) 

where F and Ba are smooth functions in the corresponding domains, Q and Sa{t,x) are a 
domain with smooth boundaries and smooth curves, respectively. Hereafter the subscripts t 
and X denote differentiation with respect to these variables, Ux'^ = = 1,2, We 

assume that BVP ([T|) and has a classical solution (in a usual sense). 
Consider the infinitesimal generator 

^ = + ^\t^ ^)^ + Vit, X, u)^, (3) 

(hereafter and rj are known smooth functions), which defines a Lie symmetry acting on 

both {t, X, u)-space as well as on its projection to {t, x)-space. Let X^^^ be the fcth-prolongation 
of the generator X calculated by the well-known prolongation formulae (see, e.g. [HE]). 

Definition 1 The Lie symmetry X of the form ^ is admitted by the boundary value 
problem (QP and ^ if and only if: 
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(a) X^^^ (^F (^x,u,Ux, . . . ,u^x''^ — = when u satisfies (Qlj; 
(h) X(sa(t, x)) = when Sait, a;) = 0, a = 1,2, . . . ,p; 

(c) X^''-'^^ [Ba (t,x,u,Ux, . . .,ui^~^^^^ = when Ba\s,{t,x)=o = 0, a = 1,2, . . . ,p. 

The definition can straightforwardly be extended on BVPs for a system of PDEs. However, 
one easily notes that this definition cannot be applied to BVPs with free boundaries, because 
such problems contain moving surfaces, say Sb(t,x) = 0, b = l,...,q, where Sb(t,x) are 
unknown functions. Obviously, these functions should be interpreted as additional variables. 
In [36] (see Appendix 2), a criteria of invariance for EVP with a free boundary was formulated. 
Another deficiency of Definition 1 appears if one consider BVPs in the unbounded domain Q 
when the boundary conditions for x = oc arise. In fact, item (b) has no sense in this case 
and cannot be replaced by the natural passage to the limit, i.e., x = L, L — )■ oo. Probably 
this deficiency for the first time was noted in [30] (see Section 4.3) where the transformation 
X = 1/y was suggested to avoid the non-regular manifold generated by x = oo. 

Now we present a definition which takes into account all possible boundary conditions and 
is applicable to a wide range of BVPs. Consider a BVP for a system of n evolution equations 
(n > 2) with 2 independent (t, x) and n dependent u = {ui, U2, ■ ■ ■ , Un) variables. Let us assume 
that the /cth-order {k > 2) basic equations of evolution type 

ui = F' {x, u,u^,..., u^^^) , i = l,...,n (4) 

are defined on a domain C with smooth boundaries. Consider three types of boundary 
and initial conditions, which can arise in applications: 

Sa{t,x) = 0: Bi {t,x,u,u^, . . . ,4''°^) = 0, a = 1, . . . ,p, j = 1, . . . ,71^, (5) 

Sbit, x) = : B[(^,x,u,. . . , Ux'\ Sb, =0, b=l,...,q,l = l,...,nb, (6) 

and 

7c(t, x) = oo : (t, X, u,u^,..., ui^"^) = 0, c = 1, . . . , r, m = 1, . . . , n^. (7) 

Here < k, k\ < k and k^ < k are the given numbers, Sait,x) and 7c(t, x) are the known 
functions, while the functions Sbit,x) defining free boundary surfaces must be found. We 
assume that all functions arising in (jl])-([7]) are sufficiently smooth so that a classical solution 
exists for this BVP. 

Consider an A^-parameter (local) Lie group Gjv of point transformations of variables [t, x, u) 
in the Euclidean space M""*"^ (open subset of R*^"*"^), which is given by equations 

t* = T(t,x,e), X* = X{t,x,e), u* = Ui{t,x,u,e), i = 1, . . . ,n, (8) 

where e = {ei, 62, ■ ■ ■ , e^) are the group parameters. According to the general Lie group theory, 
one may construct the corresponding A^-dimensional Lie algebra Ljv with the basic generators 



where 



_ dT{t,x,e) 



dsa 



(■I _ dX{t,x,e) 



£=0 



e=0 



i _ dUi{t,x,u,e) 



£=0 
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Consider the Lie algebra in the extended space ]R"+'?+2 of tlie variables {t, x, u, S), where 
S = {Si,...,Sq) are new dependent variables with respect to t and x. In the extended space 
]^n+(j+2^ the Lie group Gn corresponding to this algebra is given by transformations 

t* = T{t, X, e),x* = X{t, X, e),u* = Ui{t, x, u, s), = Sbit, x), i = 1, . . . ,n,b = 1, . . . ,q. (10) 

Now we propose a new definition, which is based on the standard definition of differential 
equation invariance as an invariant manifold Ai in the relevant space of variables and on the 
prolongation theory [3]. 

Definition 2 A boundary value problem (^-(^ is called to be invariant with respect to the 
Lie group Gn ( fT^) if: 

(a) the manifold determined by Eqs. in the space of variables (t, x,u, . . . , u*^'^^) is invariant 
with respect to the kth-order prolongation of the group G^; 

(b) each manifold determined by conditions ^ with any fixed number a is invariant 
with respect to the kath-order prolongation of the group Gn in the space of variables 
(t, x,u, . . . , u^^"^) , where ka = maxj/c^, j = 1, . . . , Ua}; 

(c) each manifold determined by conditions ^ with any fixed number b is invariant 
with respect to the kbth-order prolongation of the group Gn in the space of variables 
{t, x,u,..., M^'^'), Sb, ^) , where kb = max{kl 1 = 1,..., rib}; 

(d) each manifold determined by conditions ^ with any fixed number c is invariant 
with respect to the kcth-order prolongation of the group Gn in the space of variables 
(t, x,u, . . . , u^'^"^) , where kc = max{/c™, m = 1, . . . , ric}. 

Definition 3 The functions Ui = x), i = 1, . . . ,n and Sb = ^bit, x), h = 1, . . . ,q form an 
invariant solution of BVP corresponding to the Lie group KW\] if: 

(i) the functions Ui = and Sb = satisfy equations 0-0; 

(a) the manifold M. = {ui = x), i = 1, . . . ,n; Sb = '^b(t,x), b = l,...,g} is an 
invariant manifold of the Lie group / flOj) . 

Remark 1 Definition 2 can be generalized on more general systems (including hyperbolic and 
elliptic those) and boundary conditions containing high-order derivatives for Sb{t,x). 

Remark 2 If free boundaries are given in the form x = Sb{t), where b = l,...,q then we 
simply take Sb{t, x) = x — Sb{t) = 0. On the other hand, one can formulate a definition of Lie 
invariance for EVPs with such form of the free boundaries (see, e.g., f3Bf ). However, the form 
used in Definition 2 is more convenient for generalization on multidimensional EVPs. 

Now we present a non-trivial result to illustrate Definition 2. Let us consider the nonlinear 
BVP modeling the heat transfer in semi-infinite solid rod assuming that thermal diffusivity 
depends on temperature and the rod is exposed to a periodical flux of energy at the left 
endpoint. It should be noted that we neglect the initial distribution of the temperature in the 



4 



Table 1: Lie invariance of BVP (III])-([I3]) 



no 


k 


7 


Qo 


Lie } 


groups of invariance 


1. 


V 


V 





t* = 


te'^''^ + El, X* = xe^2+fc£3^ ^* = ue'^es 


2. 


V 





V, (go ^ 0) 


t* = 


^g(fc+2)e2 ^^^^ 3,* = xe(''+^)^^ U* = Me^2 


3. 


-2 


V, (7 7^ 0) 


V, (go 7^ 0) 


t* = 


t, X* = xe"^i, u* = we^i 



rod, i.e., consider the process on the stage when the heat transfer aheady started. Thus the 
nonhnear BVP reads as 

du d f ^du\ , , 

— = — d{u)— , t> 0, <x < +00, (11) 



dt dx \ dx ^ 
Ozt 

X = : d{u)— = gocos(7t), t > 0, (12) 

tJdj 

X = +00 : M = 0, t > 0, (13) 

where u{t, x) is an unknown temperature field, d{u) is a thermal diffusivity coefficient, go cos(7t) 
is an energy flux. We assume that all functions arising in f|TT]) -f lT3|) are sufficiently smooth, so 
that a classical solution exists for this BVP. 

Here we restrict ourselves to the case when the thermal diffusivity coefficient depends on 
the temperature as a power low, i.e. d{u) = u^, where G M, /c 7^ (in the case k = 0, the 
problem is liner and can be solved by classical methods, see, e.g., |1I]). Notably, equation ( ITTj) 
with d{u) = presents the most interesting cases of Lie symmetry invariance [Ij. In the case 
d{u) = u'^ (k G M., k ^ ~|)' admits a four- dimensional Lie group. The corresponding algebra 
L4 possesses the basic operators {dt, d^, 2tdt + xd^, kxd^ + ^udu). These operators generate the 
one-parameter Lie groups 

(14) 
(15) 
(16) 
(17) 

respectively (hereafter ei,i = 1,...,4 are arbitrary group parameters). If A; = — |, then the 
additional conformal generator (x'^dx — Sxudu) occurs, which extends L4 to the five-dimensional 
Lie algebra L5. Thus, the case k = should be examined separately. 

Theorem 1 All possible Lie groups of invariance of the nonlinear BVP / f77]) -/ f73]) with d{u) = 
, k ^ for any constants go and 7 are presented in Table 1. 

Proof. On the first step of the proof we will consider BVP f[TT]) - fll3p with the constant energy 
flux q{t) = go, i.e. 7 = 0. Let us study the case of arbitrary power k 7^ — |. First of all, 
we consider the one-parameter Lie groups (|T4|) - f|T7|) generated by the basic operators of L4. 
One easily notes that BVP (fTT|) - (fT3|) is invariant with respect to the Lie group Ti and isn't 
invariant under the Lie group T2 since the boundary curve x = isn't invariant with respect 
to the transformations fllSj) . 



Ti : 


t* 


= t + ei 




T2: 


t* 


= t, X* 


= X + 62, u* = u, 


T3: 


t* 


= te^'\ 


X* = xe^^, u* = u, 


T4: 


t* 


= t, X 


= xe^'\ u* = ue^'^ 
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According to item (b) of Definition 2, the boundary condition f lT2|) is invariant with respect 

du 

dx 



to the one-parameter group T3, if the manifold M. = {x = 0, u'^^ — qo = 0} satisfies the 



conditions 

fill* 

The first equation of (fT8l) is an identity, while the second equation leads to the expression 
goe~^^ = qo, which immediately gives 

go = 0. (19) 

The invariance of (IT^ under the one-parameter group T3 is obvious. Thus, BVP (fTT]) - (IT^ is 
invariant with respect to the Lie group T3 if and only if the restriction fll9p takes place. 

Dealing in a similar way with the Lie group T4, we obtain that BVP f|TT|) - f|T3|) is invariant 
with respect to T4 only in two cases: qo ^ 0, k = —2 and go = 0, A; G M. Indeed, according 
to item (b) of Definition 2, the boundary condition f|T2|) is invariant with respect to T4, if 
conditions (1181) are satisfied on the manifold M. Hence, we arrive at the restriction 

goe('=+2).4 ^ (20) 

which immediately leads to k = —2 provided go 7^ 0, and /c G M if go = 0. The invariance of 
(|T3|) under the one-parameter Lie group T4 is evident. 

Taking into account the restrictions considered above on go and k, one concludes that BVP 
(ITT]) - fll3l) is invariant with respect to the two-parameter Lie group Ti o T4 iff = —2, go 7^ 
and with respect to the three-parameter Lie group Ti o T3 o T4 iff /c G M, go = (it is exactly 
case 1 of Table 1). 

To find other Lie groups of invariance, one needs to consider a linear combination of the 
basic operators of L4 excepting the operator 9((we remind that the BVP is invariant under the 
Lie group Ti for arbitrary go and k) 

X = 2X-M + (A2 + (A3 + kX^)x)d, + 2X^ud^, (21) 



where Aj, i = 2, . . . ,4 are arbitrary parameters and at least two of them are non-zero. If A3 = 0, 
then one arrives only at the results obtained above for the Lie group T4, if A4 = then the 
result obtained above for the Lie group T3 is recovered. If A3A4 7^ then two possibilities occur: 
A3 -|- A;A4 7^ and A3 -|- A;A4 = 0. Consider the case A3 -|- A;A4 7^ when operator generates 
the Lie group 

Ta-. f = te^^^^\ X* = xe^^-'^''^^^'^ + , , (e(^3+fcA4)e. _ ^^ ^* ^ ue^^^'\ (22) 

A3 -f- fcA4 

Clearly, the boundary condition (fT3|) is invariant with respect to T^. Boundary conditions 
(11 2p is invariant under Ta, if and only if conditions f llSp are satisfied. Now we realize that the 
first equation of (ITSi) leads to the requirement A2 = while the second equation of ( ITS]) gives 

^^g({fc+2)A4-A3)e„ ^ (23) 

Since A3A4 7^ 0, one immediately obtains ^ = k + 2 provided go 7^ and k 7^ —2. If go = 
then we immediately arrive at case 1 from Table 1. On the other hand, the Lie group Ti o Ta 
transforms into the group Ti o T4, when k = —2. Thus, we can conclude that the BVP under 
study is invariant with respect to the two-parameter Lie group Ti o if and only if 

go^O, ^ = k + 2. (24) 
A4 
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It is exactly case 2 of Table 1. The examination of the case A3 + kX4 = leads to case 2 with 
k = —1. Thus, the invariance of BVP f|TT]) -f lT3|) with A; 7^ — | is completely examined. 

Now we examine the special case /c = — |. One easily checks that the one-parameter groups 
(with k = — |) listed in cases 1-2 of Table 1 are the groups of invariance of BVP ffTT]) - f[T^ 
(with A; = — |) under the same restrictions on the constant go- 

Thus, we need to examine whether the BVP in question can be invariant with respect to a 
Lie group corresponding to any liner combination of the basic operators of L5 

X = 2X3tdt + (A2 + (A3 + kX4)x + W)d, + (2A4 - 3A5x)u9„, A5 ^ 0. (25) 

To avoid cumbersome formulae, we consider the one-parameter Lie group corresponding to the 
pure conformal operator 

T5 : t* = t, X* = — - — , u* = (1- e^x)^. (26) 
1 — e^x 

Let us study the invariance of the boundary condition f|T3l) . According to item (d) of Definition 
2, the following equalities should take place 

x*\j>/ = +oo, u*\j\f = 0, (27) 

where Af = {x = +oc, u = 0}. However, x*\^f = — ^. Thus, the contradiction is obtained and 
we conclude that BVP ffTT]) - f[T^ with k = — | isn't invariant under T5. 

In a quite similar way, one may show that the boundary condition ( IT3|) isn't invariant under 
any Lie group corresponding to operator (125|) . 

Finally, to complete the proof, we must consider the case, when the flux of energy has 
periodical form, i.e. 7 7^ 0. Obviously, go must be nonzero, otherwise we obtain the case 
examined above. Since calculations are quite similar to the case 7 = (an analog of formula 
fl20|) plays a crucial role to derive the special power k = —2), we present the result: BVP 
f|TT|) - f|T3l) with the periodic energy flux q(t) = go cos(7t) is invariant only with respect to the 
one-parameter Lie group T4 with k = —2 (case 3 from Table 1). 

The proof is now completed. ■ 

Remark 3 Theorem 1 highlights that Definition 2 is non-trivial because the power k = —2 
isn't a special one for Lie invariance of standard nonlinear heat equation 07]) with d{u) = , 
however, k = —2 is the special power if one looks for Lie invariance of BVP / f7I]) -/ f73]) . 

3 Lie invariance of a class of (1+1)— dimensional nonlin- 
ear BVPs of Stefan type 

In this section we consider a class of (H-l)-dimensional BVPs of Stefan type used to describe 
melting and evaporation of materials in the case that their surface is exposed to a powerful 
flux of energy. Such problems also arise in mathematical modeling of other processes in biology 
(tumor growth) and physics (crystal growth). The class of BVPs after some simplifications 
(like using the Goodman substitution to transform the basic equations to the standard heat 
equations) can be written as follows 



dv 
dt 



dx \ dx J ' 



(29) 



5*1 (t, x) = : di{u) 



du 
dx 



H^{u)Vi-q{t,u), Vi = h{t,u), 
du 



(30) 



S2{t,x) = : d2{Vm) 



dv 
dx 



(31) 
(32) 



X = +00 : V = V, 



where u{t, x) and v{t, x) are the unknown temperature fields; Sk{t, x), k = 1,2 are the unknown 
functions, which determine the phase division boundaries (they can be also presented in the form 
Sk{t,x) = X — Sk{t)); Vk{t,x) = —^§f/^§^ k = 1,2 are the phase division boundary velocities; 
q{t,u) is the known strictly positive function presenting the energy fiux being absorbed by 
the material; h{t, u) is the known non-negative function describing dynamics of evaporation 
process; H^, k = 1,2 are the known strictly positive function presenting specific heat values 
per unit volume of liquid and solid phases. The parameters Um, Vm and f oo are assumed to 
known, moreover, Vm Voo- 

Here Eqs. fl28|) and fl29|) describe the heat transfer process in liquid and solid phases, the 
boundary conditions fISUl) present evaporation dynamics on the surface Si{t,x) = 0, and the 
boundary conditions fl?T]) are the famous Stefan conditions on the surface S2{t,x) = dividing 
the liquid and solid phases. Assuming that the liquid phase thickness is considerably less than 
the solid phase thickness, one may use the Dirichlet condition (l32i) . It should be stressed that 
we neglect the initial distribution of the temperature in the solid phase and consider the process 
on the stage when two phases take already place. 

One may claim that formulae fl28l) - fl32l) present a class of EVPs with moving boundaries and 
take into account a number of different situations, which occur in the melting and evaporation 
processes. Setting q(t, u) = const, Vi = const$(t) and h{t, u) = ^(t)u, where $(t) is a correctly- 
specified function, one obtains the problem, which is the most typical, see, e.g., [30]. In the case 
of a process when surfaces are exposed to very powerful periodic laser pulses these functions 
take complicated forms |39] . 




The EVP obtained is based on the standard nonlinear heat equations. Lie symmetry of non- 
coupled system (I28|) - (l29|) can be easily derived using the determining equations from paper [i2] . 
where reaction-diffusion systems of more general form have been investigated. Now we formulate 
a theorem, which gives complete information on Lie symmetry of this system. 

Theorem 2 All possible maximal algebras of invariance (up to equivalent representations gen- 
erated by transformations of the form ) of the system and for any fixed vectors 
{di,d2) with strictly positive functions di{u) and d2{v) are presented in Table 2. Any other 
system of the form l[2^} and [29^) is reduced to one of those with diffusivities from Table 2 by 
an equivalence transformation of the form 



where 7^ (i = 0, . . . , 3), to, ^^o, Uq, and Vq are arbitrary parameters. 

Remark 4 In the case of linear system [2B^) - [29\) with ki = k2 (see case 9 of Table 2), the 
Lie algebra extension occurs by the operators vd^ and ud^. However, EVP / f^) -/ f5'^) with 
di{u) = d2{v) = ki is rather artificial from physical point because diffusivities of solid and 
liquid phases must be different. Thus, we don't consider this case below. 




t— J- Cot + to, X ^ eix + Xq, u^e2U + UQ, v ^ e^^v -\- Vq, 



(33) 



8 



Table 2: Lie algebras of NHEs system (l28l) - (l29|) .(Here ki,k2,m and n are arbitrary non-zero 
constants; while a{t,x) and (3{t,x) are arbitrary solutions of the linear heat equations at = 
kia^x and Pt = k2f3xx, respectively.) 



no 


di(u) 




Basic operators of MAI 




1. 


V 


V 


A= {dt,dx,2tdt + xdx) 




2. 


ki 


V 


A, udu, oi(t, x)du 




3. 


V 


k2 


A, vdv, /3(t, x)dy 




4. 


e" 




A, xdx + 2du + 29„ 




5. 


e" 




A, xdx + 2du + ^vd^ 




6. 






A, xdx + "iudu + 2(9^ 




7. 






A, xdx + ludu + |^w<9^ 




8. 


4 

M 3 


4 

V 3 


AjXC^i: — \udu — \vdv,x'^dx — Sxudu - 


- ?>xvdu 


9. 




k2 


A, udu, vd^, a{t, x)du, (3{t, x)d^, 2tdx - 
4txdx + U^dt - ^(x^ + 2kit)udu - ^ 


X i^-^udu + ^^vd,^ , 

(x^ + 2k2t)vdu 



Table 3: Lie invariance of BVP ()28|)-fl32|l 

no q{t,u) h(t,u) Lie groups of invariance 

t* = t, X* =x + e, u* = u, V* = V, SI = Si, S^ = S2 

t* =t + El, X* = X + 62, u* = U, V* = V, SI = Si, S2 = S2 

t* = te^^\ X* = xe^i + 62, u* =u, v* = v, SI = Si, S^ = S2 



L V V 

2. q{u) h{u) 

q £W h(u) 



Remark 5 // one takes into account the trivial discrete transformations t t, x x, u ^ v, 
and f — 7- u, then cases 2 and 3, 5 and 6 arising in Table 2 are equivalent. However, the class 
of EVPs [2^) - [3^] isn't invariant under these transformations because of boundary conditions 
[3D\) and [3^) . Thus, we don't take into account them in what follows. 

Using the set of equivalence transformations fl33l) . we can straightforwardly extend one to 
the relevant set for BVP fl28|) - fl32|) by adding the identical transformations for the variables 
Sk{t, x). Direct calculations show that the most general form of those is 

t^eot + to, X ^ eix + xo, u e2U + uo, v ^ e^v + vo. Si ^ Si, S2 ^ S2. (34) 

where Cj 7^ {i = 0, . . . ,3),to, xq, uq, and vq are arbitrary parameters (ei > 0). 
Now we formulate the main result of this section. 

Theorem 3 BVP [2^) - ^3^) with any smooth functions di{u), d2{v),q{t,u),h{t,u) and Hi{u) 
is invariant under the one-parameter Lie group presented in case 1 of Table 3. All possible 
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extensions of this Lie group invariance (up to equivalent representations generated by equiva- 
lence transformations of the form [34\)) depend only on the form of the functions q{t,u) and 
h{t,u), and are presented in cases 2 and 3 of Table 3. Any other BVP of the form ^2^) - [3^) 
is invariant under two-parameter Lie group is reduced by transformations (3^ to one of those 
with the functions q{t, u) and h(t, u) from Table 3. 

Proof. According to Definition 2 and Theorem 2 we need to examine the nine different cases 
hsted in Table 2. It turns out that the examination of the first case, when the functions di{u) 
and d2{v) are arbitrary, leads to the main result of the theorem presented in Table 3. 

Let us consider the one-parameter Lie groups corresponding to the basic operators of algebra 
A. Obviously, BVP (I28l) - (l32|) with arbitrary given functions is invariant under the group of 
space translations generated by the operator = dx and this is listed in the first case of 
Table 3. Since any linear combination XiPt + A2-D of other two operators is equivalent (up to 
transformations either to Pt = dt (if A2 = 0) or to D = 2tdt + xd^ (if A2 7^ 0), we should 
separately examine these two operators. 

Now we apply Definition 2 to D . Taking into account that BVP fl28|) - fl32|) has two free 
boundaries, we construct the extended Lie group To corresponding to the operator D: 

fo-.t* = te^'\ X* = xe'\ u* = u, v* = v, SI = S; = S2. (35) 

According to item (c), the boundary conditions fl30|) are invariant with respect to the group T^, 
if the manifold M = {Si{t,x) = 0, di{u)^ = Hi{u)Vi - q{t,u), Vi - h{t,u) = O} satisfies the 
conditions 



Sl\j^ = 0, d^{u*)— - H^iu*)V* + q{t\u*] 



du* 

Taking into account fl35|) . one finds 



0, V* -h{t\u*)\j^ = Q. (36) 



M 

M 



dx* ^ dx'' dx* ^ dx' ^ ^? 55 (37) 

so that the second and third equations of fl5B]) produce the equations 

e''q{te^'\u) = q{t,u) e''h {te^'\u) = h{t,u), (38) 

to find the functions q{t,u) and h{t,u). Solving (138|) one obtains 

qM , , , h(u) 
q{t,u) = ^, h{t,u) = ^, (39) 

where q{u) and h{u) are arbitrary smooth functions. The invariance criterium of the boundary 
conditions fl3ip for is fulfilled for arbitrary parameters arising in fl3ip . while the invariance 
of condition ( 13^ under Td is obvious. Thus, BVP is invariant with respect to the 

Lie group if and only if restrictions ( l39|) take place. This is exactly listed in case 3 of Table 
3. 

In a quite similar way one can show, that the BVP under study is invariant with respect to 
the extended Lie group corresponding to the operator Pf = dt if and only if the restrictions on 

q(t,u) = q{u) and h{t,u) = h{u), (40) 

take place, and this is what exactly listed in case 2 of Table 3. 
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Much more cumbersome calculations are needed to show that there are no any new Lie 
group invariance for BVP ( I28|) - fl32|) nevertheless there are eight special cases listed in Table 2, 
which lead to the extensions of MAI of the basic equations fl28|l . 

Let us consider case 2 of Table 2. Firstly, we check the invariance of BVP f l2S]) - fl52]) with 
respect to the one-parameter extended Lie groups corresponding to the operators Xi = udu 
and X2 = a(t,x)du: 

fi-.t* = t, X* = x, u* = ue'\ V* = V, SI = Si, S; = S2, (41) 

and ^ 

T2 : t* = t, X* = X, u* = u + a(t, x)e2, v* = v, SI = Si, S2 = 5*2. (42) 

According to item (c) of Definition 2, the boundary conditions fl3T|) are invariant with respect 
to the group Ti, if the conditions 



= 0, U* - Um\j^ = 0, V* - Vm\j^ = (43) 

AT 



dv* du* 
are satisfied, where the manifold 

M = < S2{t,x) = 0, d2{v„,)— - di{Um)-^ - H2{Vm)V2 = 0, M - M„ = 0, V - = > . 

Taking into account fj^Tjl and the second equation of fj43|) . we arrive at the requirement 

du du ^ , , 

— = —e'^ ^ei = 0. 44 
ox ox 

Similarly, one easily checks that the boundary conditions (I3T!) isn't invariant with respect to 
the Lie group T2, too. Indeed, to satisfy the third equation of fH51) . one obtains the requirement 

a{t, x)e2 = ^ £2 = 0. (45) 

Let us now examine the invariance of BVP with respect to an extended Lie group 

Tc corresponding to a liner combination of operators Pt, D, Xi, and X2, i.e. 

X, = (Ai + 2\2t)dt + AsxS,. + (A3M + A4a(t, x))du, (46) 

where Aj, i = 1,...,4 are arbitrary parameters and A| + A4 7^ (otherwise the operator 
XiPt + A2-D is obtained). Having transformations (IMl) . we can put Ai = and A2 = 1 in (B6l) 
so that the operator takes the form 

X, = 2tdt + xda, + (A3M + \ia{t, x))du, A^ + A^ ^ 0. (47) 

The corresponding Lie group T^. is 

f,:f = X* = xe'% u* = ue^''^ + A4 /" a{te^\xe^)e^'^^-'^^dT, v* 

Jo 

SI = Si, S*2 = S2. (48) 

Now we again show that boundary conditions (ISTj) are not invariant under T^. In fact, the third 
equation of fH5]) gives the restriction 



Sc 



A4 / a{te^\ xe")e^^("-"^)c/r = m„(1 - e^^"^, (49) 
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so that the Lie group Tc can be written in the form 

fe : r = X* = u* = ue^''^ + - e^^"^), v* = v, S*^ = Si, S; = S2. (50) 

Taking into account fl50|) and the second equation of fj43|) . we arrive at the requirement 

= ffe^- ^ A3 = 0, (51) 

OX ox 

what leads to A4 = (see fHU]) ). However, this contradicts to the assumption A3 + A4 7^ 0. 

Hence, case 2 of Table 2 is completely examined. Cases 3 and 9 of Table 2 can be studied 
in a quite similar manner because each of them needs to examine groups (HT]) and (142|) . 

Consider case 4 from Table 2. Here the operator X4 = xdx + 2du + 2dy arises, which generates 
the extended Lie group T4: 

t* =t, X* = xe^\ u* = u + 2^4, V* = v + 2^4, = Si, S*^ = S2. (52) 

Applying Definition 2 to the boundary conditions (I5T]) in the case of T4, equations (HSj) are 
again obtained. The third and fourth equations of ( H3|) lead to the requirement 

Um=Um + 254, Vm = Vm + 2^4, (53) 



hence, £4 = 0. Thus, EVP fl28|) -f l32|) cannot be invariant with respect to the Lie group T4. 
Moreover, the same result is obtained if one examines any linear combination of operators Pt, 
D, and X4, i.e. 

X = (Ai + 2\2t)dt + (A2 + \3)xd, + 2\3udu + 2\svd,, A3 ^ (54) 

with A3 7^ 0. Thus, we conclude, that the exponential diffusivities arising in case 4 of Table 2 
don't lead to any new Lie groups of invariance of EVP ( I28|) -( l32|) . Cases 5 and 6 of Table 2 can 
be studied in a quite similar way. 

Let us consider case 7 of Table 2, which needs separate examination. The extended Lie 
group T7 corresponding to the operator X7 = xd^ + ^udu + ^vd^ is 

Tj:t* = t, X* = xe'\ u* = uei'\ v* = ve^'\ S^ = Si, S; = S2. (55) 



In order to the boundary conditions fl3Tl) be invariant under T7 equations fl43|l are again 
obtained. It turns out, equations fH3|) are fulfilled ifm>0,?T,>0, = Vm = and -^2(0) = 0. 
Hence, we must apply item (d) of Definition 2 to the boundary conditions 



x*\^ = +oo, v*-v^\^ = 0, (56) 

where V = {x = +00, v — v^o = 0}. The group Tj transforms the second equation of fl56l) as 
follows 

VooC^'' -v^ = 0, (57) 

hence, we arrive at the restriction foo = (otherwise £7 = 0). Thus, the contradiction is 
obtained because, in the very beginning, we assumed Vm 7^ v^o- One may check that the same 
result is obtained for any liner combination of the operators Pt, D, and X7 . Thus, case 7 of 
Table 2 is completely examined. 

Case 8 of Table 2 can be treated in a quite similar way as we did in Theorem 1 (see formulae 
(126]) and dSZD). 

The proof is now completed. ■ 

Finally, we note that theorem 2 from the recent paper [57] follows as a particular case from 
Theorem 3 (but not wise versa !), nevertheless Definition 2 was not used in 
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4 Symmetry reduction and invariant solutions of the 
class EVPs ([28D-([3D 



In this section we consider symmetry reduction of BVPs of the form f l25]) - fl52]) to BVPs 
for systems of two ordinary differential equations and construct exact solutions for the reduced 
BVPs. 

Operator P^, corresponding to the invariance transformations in case 1 of Table 3, leads to 
an ansatz, which doesn't depend on the space variable x. This contradicts to the free boundary 
surfaces and leads to non-physical solutions. 

According to case 2 of Theorem 3, each BVPs belonging to the class under study for q{t, u) = 
q{u) and h{t,u) = h{u) admits the two-dimensional Lie algebra with basic operators Pt = dt 
and Px = dx- Hence, it also admits the operator Xi = dt + fidx, /i G M, which leads to the 
plane-wave ansatz 

u = V = Sk = SkiOi i = x- lit, (58) 

where k = 1,2 and fi is an unknown velocity. Note, that ansatz fl58l) with /i = leads to 
stationary solutions of the BVP in question. However, these solutions don't have essential 
physical sense so that will not studied hereafter. 

In case 3 of Table 3, each BVP of the form fl28|) - fl32|) admits two-dimensional Lie algebra 
with basic operators Px and D = 2tdt + xdx- Obviously, that any linear combination of these 
operators is equivalent, up to transformations flM|) . to operator D, which generates the ansatz 

X 

u = u{uj), V = v{uj), Sk = Sk{oo), 00 = k = 1,2. (59) 

Using ansatze (lS5]) - flSU]) and taking into account the restrictions on the functions q{t, u) and 
h{t, u) (see Theorem 3) one can easy obtain the following theorems. 

Theorem 4 Ansatz / f5^) reduces any nonlinear BVP of the form [2B^) - [3^] with the coefficient 
restrictions q{t,u) = q{u) and h{t,u) = h{u) to the BVP for the second-order ODEs 



d f , , . du\ du ^ ^ , , 

d f , , ,dv\ dv ^ ^ ^ X 
dzL 

e = : di{u)— = Hi{u)fi - q{u), /i = h{u), (62) 

^ = 5 : d2{v.m)^ = dl{Um)^ + i^2(^^m)/i, U = Um, V = V^, (63) 

^ = +00 : V = Vq, (64) 
where 6 and fi are to be determined parameters. 



Theorem 5 Ansatz (E^) reduces any nonlinear BVP of the form ^2^-^3^ with the coefficient 
restrictions q{t,u) = and h{t,u) = to the BVP for the second-order ODEs 

d f , , .du\ uj du 

d f , , .dv\ UJ dv 
d^i^^^"^d^J + 2rf^ = °' 
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Ui < (jJ < UJ2, 



(65) 
(66) 



uj = ui: di{u)'^ = Hi{u)y - Qiu), y = h{u), (67) 
, , .dv du UJ2 

UJ = UJ2. "2(fm)^ = "H""^)^ + H2{Vm)Y' ^ = ^m, V = Vm, (68) 

oj = +00 : V = Vq, (69) 

where Ui and 002 are to be determined parameters. 

Let us consider BVP fl60l) - fl64l) . It is well-known, that the system of nonlinear ordinary 
differential equations fl60|) and fl6T]) with fi = can be linearized by the Kirchhoff substitution. It 
turns out, this substitution can be generalized in the case /i 7^ by introducing new independent 
variables |30j 

U = U- Um, V = V -vo, (70) 

^= r d,{U + Um)dr], 0<^<(5, (71) 
^0 

C = 6 + £d2{V + vo)dr], C>S, (72) 

where the lower bonds of integration are chosen as ^|^=o ~ ^ ^\'q=s* ~ ^■ 

Substituting ([70D-([72]) into BVP ([60D-([64l), one obtains the BVP for system of Unear dif- 
ferential equations 

d'^U dU 

_ + ^_ = 0. 0<„<S', (73) 
dV „ 

_+^_ = 0, „>^-, (74) 

= : ^ = Hi{U + Urr,)lJ^ - q{U + u„), U = h'^{^) - Um, (75) 
drj 

V = S*: ^ = ^ + H2{vm)^, U = 0,V = Vm, (76) 

ri = +00: V = 0, (77) 

where is an inverse function to h{u) (hereafter existence of this function is assumed), 
Vm = Vm — Vq, 6* is ucw to bc determined parameter. 

To solve BVP ( 175]) - (177ll one needs to find the unknown functions U = Uiji) and V = V{r]) 
and the parameters fi and S*. Since equations f l73|) and ([71]) are linear ordinary differential 
equations with constant coefficient their general solution is 

U = Ci + ^26"^", 0<r]<5*, (78) 
^ = C3 + C4e-^^ r7>5*, (79) 

where Ci{i = 1, . . . , 4) are to be determined constants. 

Substituting solution (178|) into the second equation of the boundary conditions ( I75|) and the 
second equation of ( 176|) . one finds the constants Ci and C2: 

Ci = {h-\fi) - Um) ^^J!^, C2 = {h-\fi) - Um) ^^^y (80) 
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Similarly, the constants C3 and C4 can be found using the third equations from fj76|l and the 
boundary condition fl77|) 

C3 = 0, C^ = Vme^'\ (81) 

Thus, using formulae (178|) -( 18T|) . we construct the explicit formulae for the function U = Uirf) 
and V = V{ri): 

fi{S*-ri) _ 1 

U = {h-\f,) - u^) , (82) 

V = Kne'^(^*-''\ (83) 

Finally, we need to specify the parameters fi and 6*. This can be done by substituting (15^ and 
fl83|) into the first equations of the boundary conditions f l75l) and fl76|l . and using the formulae 

^^-.r'M-...)^^;?^. ^ - (84) 

After the corresponding calculations, we arrive at the equations 

{h-\l^) - um) = ^^^^^ - H,ih-\^)), (85) 

Vm = (h'^ifl) - Um) ^f,s*^_ ^ - H2{v.m)- (86) 

Thus, equation 0861) leads to the explicit formula for the parameter 5*, which corresponds to 
the thickness of liquid phase: 

/i Vm + H2iVm)J 

Having fl87|l . equations fl85l) produces immediately the transcendent equation to find the velocity 
^" _i 

^a^^M _ H,{h-\^)) - h-\fi) = Vm-Um + H2{Vm). (88) 

/i 

Thus, formulae ( I82|) - (l83|) give exact solution of problem ( [73|) -(I77I). where the important 
parameters S* and /i are defined by expressions (!87|) and (!88l) . respectively. It means that 
the exact solution of BVP fl60p - fl64p is obtained in implicit form. It should be stressed that 
this result is essential generalization of [20], where a particular case of BVP was 
investigated. 

Example. We use a model, which was formulated and investigated by numerical methods 
in [3H] • The model describes the processes of heating, melting and evaporation of metals under 
the action at their surface of powerful laser pulses. Under the relevant assumptions it can be 
written as follows 



d ( dT \ dT 

^ f Ai(Ti)^j =pci(Ti)^, < si(t) < X < S2(t), t > 0, (89) 

5j(A,ffi)^)=pc.(T,)-^. .>.,W.<>0. (90) 

dT It / T \ 
x = s,{t): Ai(Ti)^ = pi:,ri-Q(t, TO, V^i = Ky^ expf-^j, (91) 

BT BT 

X = S2{t) : HTm)^ = + pLmV2, T, = T^ = T^, (92) 

X = +00 : T2 = Too, (93) 
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where T^, T^, Too are the known temperatures of evaporation (under normal atmospheric 
pressure), melting and solid phase of metal, respectively; Afc(Tfc), Cfc(Tfc), p, and are 
the specific heat coefficients (functions), which are typical for the given metal (note, that we 
consider the model with constant and equal densities of solid and liquid phases of metal, i.e., 
pi = p2 = p); Sk{t) are the phase division boundary coordinates to be found; Vkit.x) 



IT 



are the phase division boundary velocities; Tk{t,x) are the unknown temperature fields; index 
k = 1,2 corresponds to the liquid and solid phases, respectively. 

We will consider the processes of melting and evaporation under a long time powerful pulse, 
i.e. ~ 10~^ sec and more. In this case, the function Q{t,Ti) determining the amount of the 
absorbed energy is defined by the formula 

Q(t,Ti)=x(Ti)-go(t), (94) 

where qo(t) is the power of laser pulse, assumed to be a constant go, and x(Ti) is the absorbtion 
coefficient of the energy. Velocity which is approximately equal velocity of sound in the 
metal, can be given by the formula [39] 

P„yi T. AL, 

where Pa is the normal atmospheric pressure, A is the atomic weight, R is the universal gas 
constant. 

Since our aim is to demonstrate that the formulae obtained above produce realistic data, we 
need to specify all coefficients arising in BVP f l5^ - fl^?5]l . Thus, all coefficients were taken from 
the paper [32], where the processes of melting and evaporation of aluminium were studied. In 
the particular case, the temperature dependence of the quantities CkiTk) and x(Ti) are essential 
and have the form 

Ci(Ti) = 1086, T > T^, C2(T2) = 752.2 + 0.473 ■ T, T < T^, (96) 

- rp \ 0.4 
1 



where [cfc(Tfc)] = Jkg^^K^^. Other physical quantities are assumed to be some constants, 
namely: Ai = A2 = 240 Wm-iK"\ p = 2545 kgm'^, = 0.64-106 Jkg-\ L„ = lO.S-lO^ Jkg-\ 
= 2793 K, Tm = 933 K, = 300 K. 

First of all, to use the results obtained above, we should transform the governing equations 
of BVP dMD-dnSD to the form and Using Goodman's substitution 

Ti T2 

u = 0i(Ti) = j ciOpdC, V = 02(T2) = J c^iOpd^- (98) 


one easy transforms BVP (189D-([93D to the form (128])-([32D, where 

d,{u) = ^, d,{v) = ^ ■ , (99) 

0-64 go ,0.4 PaVciA __T^ 1 pciT^ 



^(") = ii6oo^-(;^"'' ^(") = 75^"^^'7^''H"^'' ^'''^ 
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and 

H, = pL,, H2 = pLm. (102) 

In formulae ( I99|) - fll02p . we use the notation C2(T2) = 02 + ^2^2, where the coefficients 02 and 
62 are determined by (196|) . 

Now, using formulae (182|) and (!83|) . and (l87|) and (!88|) . and making the relevant simplifica- 
tions, one can receive the exact solution of BVP f l89p - (IU5]) in the explicit form 

pCi 

y ^ T00C2 (^) + (T^ - T^)C2 (^) e-"^^^'^) 
C2(^)-|(T.-T.)e-^(^-) 

5=^logfl + ^^:^^^^^^^V (105) 

where K = {Tm — Too)pc2 ( '^'"^^°° ) + if2 and the velocity p satisfies the transcendent equation 

q{h-\p)) 



P 



h-\p)=K - pciT^ + Hi. (106) 



Equation fll06p can be easy solved numerically by means of Maple (Mathematica etc.) 
program package. We used Maple 12. The calculations were carried out for two values of the 
parameter go- 1) Q'o = 1 ■ 10^° Wm~^; 2) go = 5 ■ 10^*^ Wm~^. The following phase division 
boundary velocities {[p] = m • c~^) were obtained: 
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0.10 if go = 1 ■ 10 
^ ~ 0.54 if go = 5 ■ 10^°. 



The corespondent liquid phase thickness {[5\ = m) is 

_ e _ / 9.60 ■ 10-^ if go = 1 ■ 10^°, 
S2[t) si[t) - - I 2.23 ■ 10-4 if go = 5 ■ lO^". 

The temperature fields of liquid and solid phases of aluminium are presented in Fig.l. Compar- 
ing the temperature fields and the liquid phase thickness obtained here with those from [39], 
one concludes that they are sufficiently similar. Of course, one should take into account that a 
chain of laser pulses was used in for numerical simulations, while one only laser pulse was 
used to obtain formulae presented above. 



5 Conclusions 

In this paper, nonlinear boundary value problems by means of the classical Lie symmetry 
method are studied. First of all, an analysis of the known definitions of Lie invariance for 
BVPs is presented. Having this done, a new definition of invariance in Lie sense for BVPs 
of the form (ll])-([7]) is formulated. This definition (see Definition 2) is applicable to a wide 
class of BVPs, including those with several basic equations, with moving boundaries, and with 
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Figure 1: Exact solutions of problem f l89|) -f l93|) for aluminium with the energy flux: 1) go = 
1.0 ■ 10^° Wm-2 and 2) go = 5.0 ■ 10^° Wm'^. 



boundary conditions on non-regular manifolds, moreover, it can be easily generalized on BVPs 
with hyperbolic and elliptic basic equations. 

The class of two-dimensional nonlinear BVPs ( I28p -( 15^ . modeling the process of melting 
and evaporation of metals (under acting a powerful flux of energy) is studied in details. Using 
Definition 2, all possible Lie symmetries (see Theorem 3) and the relevant reductions with 
physical meaning to BVPs for ordinary differential equations (see Theorems 4 and 5) are con- 
structed. The example how to construct exact solution of the nonlinear problem fl25]) -f l5^ with 
correctly-specified coefficients for aluminium is also presented. We established that exact for- 
mulae fll03p -( IT05i) . obtained by direct application of Theorem 4, lead to quite realistic results, 
which are sufficiently similar to those obtained earlier by numerical simulations. 

The work is in progress to extend the results on multidimensional BVPs using the definition 
proposed in this paper. 
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